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ABSTRACT 


In this paper the transient thermal stress problem in a hollow cylin- 
der or a disk containing a radial crack is considered. It is assumed 
that the cylinder is reinforced on its inner boundary by a membrane which 
has thermoelastic constants different than those of the base material. 

The transient temperature, thermal stresses and the crack tip stress 
intensity factors are calculated in a cylinder which is subjected to a 
sudden change of temperature on the inside surface. The results are 
obtained for various dimensionless parameters and material constants. 

The special cases of the crack terminating at the cylinder-membrane inter- 
face and of the broken membrane are separately considered and some exam- 
ples are given. 


1 . Introduction 

Cracking due to thermal stresses arises in many practical applica- 
tions particularly when sudden changes occur in the environmental temper- 
ature. Suddenly cooled glass plates and containers are generally the 
familiar examples for the phenomenon. The problem may also be quite impor- 
tant in such structural components as the pressure vessels, piping and 
hollow circular disks subjected to thermal transients. In these structures 
the actual problems are usually very complicated three-dimensional crack 
problems. Under thermal shock a part-through crack may initiate or an 
existing pa rt-th rough crack may propagate in axial or circumferential as 
well as in the thickness direction. However, if the crack driving force 
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due mechanical stresses is sufficiently small, then because of the self- 
equilibrating nature of the thermal stresses the growth of the crack in 
thickness direction is generally arrested. A clear demonstration of this 
process may be seen in the results reported in [1] where it was observed 
that in relatively thick-walled hollow circular glass cylinders suddenly 
cooled from inside an axial initial crack penetrated into the cylinder 
wall only partially but propagated axially along the entire length of 
the cylinder. It is clear that in such problems the corresponding plane 
strain or plane stress solutions may provide very useful bounds for the 
actual part-through crack problem. 

The axi symmetric problem for a thick-walled cylinder containing a 
part-through circumferential crack and subjected to transient thermal 
stresses was considered in [2]. In this paper we consider the correspond- 
ing plane strain or plane stress problem for a hollow circular cylinder 
which contains a radial crack. To simulate the cladding (in, for example, 
nuclear pressure vessels) it is assumed that the cylinder is reinforced 
on the inside boundary by a (thin) membrane. The results found may also 
be applicable to other composite cylinders or disks in which the thick- 
ness of the inner cylinder or ring is relatively very small. 

In solving the problem it is assumed that the material is linear, 
the thermo-mechanical constants of the cylinder and the reinforcing mem- 
brane are independent of temperature, the bending stiffness of the rein- 
forcing membrane and all thermoelastic coupling effects are negligible, 
and the transient thermal stress problem may be treated as a quasi -static 
time-dependent problem, that is, all inertia effects may be neglected. 
Previous studies on dynamic thermoelasticity indicate that this last 
assumption, which simplifies the analysis of the problem quite consider- 
ably, would not cause any significant changes in the results (see, for 
example, [3] and [4]). Since the problem is linear the solution of the 
crack problem due to thermal stresses and to all other sources of load- 
ing may be considered separately. Also, the solution of the crack problem 
under thermal stresses may be obtained in tv/o steps. The first would 
be the calculation of time-dependent thermal stresses in an uncracked 
reinforced cylinder under given temperature and stress boundary conditions. 
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Then the second step would be the solution of an isothermal crack prob- 
lem in which the self-equilibrating crack surface tractions equal and 
opposite to that obtained in the first solution are the only external 
loads. The complete solution may be obtained by adding these 1 uS t two 
solutions to the other isothermal results. Needless to say, in thermal 
stress as well as in the isothermal problems the important information 
which is useful in fracture considerations is contained in the perturba- 
tion solutions. 


2. The Thermal Stresses 

Referring to Figure 1 for notation let the temperatures in reinforcing 
and main cylinders be T-j and respectively. If we assume that the 
crack surfaces are well -insulated, then T-j and I2 would be functions of 
radial coordinate r and the time t only. Let Tq be the initial tempera- 
ture of the composite cylinder and the temperature changes e-j and Sg be 
defined by 


e^(r;t) 


, (a^j<r<a, O^t) , 

(1) 

e2(r,t) 

= TjCr.t) - 

s (a<r<b, 0<t) . 

(2) 


Ir the case of uncracked cylinders the problem is axisymmetric, the shear 
stresses are zero, and the relevant displacement and stress components 
may be expressed as follows (see, for example, [5]): 

9i(p.t)pdp + C^.r + -^ , (3) 

’ a. 

1 



“ 1^1 


1 


E_. 


e^(p,t)pdp T (1-2L 


®i 



(4) 


aJi(r,t) 


^ 1 
1-v^. 


E-a • 


E. 


e,(p.t)pdp -^e,(r.t) 

1 11 

(5) 
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li 


where i = 1,2, 


, for i = 1 , 

a^. = a , a<r<b , for i * 2 , 

a., E. and are the standard thernioelastic constants of the materials 
is the radial displacement and the unknown constants and C2i *'•« 
determined from the following boundary and Cv< ! tinuity conditions: 


aJl(ao,t) =0 , ® > 

u|^(a,t) = Ug^lajt) , aj-|(a,t) ~ * 

Thus, from (3), (4), (6) and (7) it can be shown that 
1 -2Vn Am A-i 



a«(l+Vp)(l-'2v«)B 


+ 


1-2v2 A2 

-w-i; 



(6a ,b) 
(7a,b) 


(8a-d) 


E2(a2-b^)[(1-2vi)a2+a2] E,(a2-a^^)[(1-2v2)a ^+b^] 

*0 ” ‘ n+vJJa^^r*F^ * (l+vi)a^^a‘*b^ ’ 

E,(a2-b2)[-(l+u,){l-v.)a,b’A + {1-u,)(l+U2){l-2u,)a2a2B] 
'^l (l-v,)(1-v2‘!)i'‘Fr- 

[a2E2a^(l-v^ )B + aiE^b2(l-«2)A][(l-2v2)a^+b^] 

+ n-v,)(1-V2)a‘*b^ 


9 
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[(l-2v^)a2+aQ2][j4^^^a2(T.^^)B + a^E^b2(l-v2)A] 

(i-v,)(l-v2)ao'^a‘'b2 


+ 


(a2-a^^)[(1+v^)(l-V2)a^b2A - 0-v^){l+V2)(l-2v2)a2®^^3 
rr^TTrv-^T^-^b^ ^ 


t 



(9a-e) 


It should be noted that the expressions (3)-(5), (8) and (9) are valid 
for the case of plane strain only. For "disk" problems which may be 
approximated by generalized i.lane stress E. and should be replaced by 
E^.(l+ 2 v^.)/(l+v^ 2 ) and v^/(l+v^. )» respectively. 

To obtain the temperature distributions 9^ and 02 the following dif- 
fusion equations must be solved under given initial and boundary condi- 
tions; 

i 894(r,t) 

v2e,(r,t) = L , (1=1,2) (10) 

where aQ<r<a for i * 1, a<r<b for i = 2 and and D 2 are the respective 
coefficients of diffusivity. A particular problem of interest is the 
sudden cooling of the composite cylinder from inside which is initially 
under a homogeneous temperature T^. Thus, if it is assumed that for t>0 
the cylinder wall is maintained at a constant temperature T^^ and 

the outer surface r = b is insulated, (10) must be solved under 


8^(r,0) = 0 , (aQ<r<a) , (11) 

02(«"»O) == 0 , (a<r<b) , (12) 

0^(a,t) = 62(3, t) , -^.“©^(ajt) = (I3a,b) 

0^(ag,t) = 0aoH(t) , ^02^^’*^ * ° (14a,b) 
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where 
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ao 
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(15) 


and kg are the coefficients of heat conduction and H(t) is the 
Heaviside function. By using the Laplace transforms the solution of 
(10) subject to (11)-(14) may be obtained as follows: 


ei(r*.t*) „ 2 


-X 2t* 


a*0*(x„) ^S^^oo^®*^n * ^*^n ^^1 1 ^ * ^^n ^ 


+ Zo,(6a*x„.SX„)2^,(r*x„,a*x„)] , 

92(r*,t*) _ 4 I " 

®ao ” 1 


(16) 

(17) 


where the cylinder functions and the dimensionless quantities are defined 
by 


Z^j(x,y) = J^.(x)Yj(y) - Y.(x)Jj(y) , (i,j=0,l) , (18) 

r* = r/b , a^,* = a^/b , a* = a/b , t* = D,t/b2 , s = k^, = k 2 /k■^ . 


(19) 

Z„,(6X„.6a\)] 

- SKo^n(8Xn.Ba\)[-Zjj,(a|j*x„,a*x^) + "^ Zjj^(a*Xj,,a|j*Xj,)J 

* S^ol(®o*^n’^*^n>C-^ll<s^n’^"*^n> ^ . (ZO) 

are the roots of 
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( 21 ) 

and and are the Bessel functions of first and second kind, respec- 
tively. 

The thermal stresses in the composite cylinder without the crack 
may be obtained by substituting from (16) and (17) into (4) and (5). 

For example, the hoop stresses and which are needed in the solu- 
tion of the crack problem are found to be 


-y. - ■ C-ptz 


aO 




(22) 


,T , . ... . “2^2830 / 2 <'’*) <i,(l+r* 2 ) ^ __ 

’02''^ (“Fsz eT:; ^2*’^ '1 , -r*’^ — , 


where 


I 

(r*) = I 


e^(r*,t*) 


ao 


ao 


.. . 

,*dr* = —2 . 


(23) 


( 6a*x„ . ex „ ) Cr*Z„, ( an„ , r*X„) -a^*Zo, (a*^„ .a„n„ ) ] 
+ Z„,(sa*A„,6X„)[r*Z,i(rn„,an„)-a„*Z„(aJx^,a*X„)]} , 



02(i"*.t*) 


ao 


r*dr* 


-\ ^t* 

- r*^-a*^ 4. 4 ” e 

2 TT 5 3a*^X^^H*(x^T 


•[r*Zp(er*X^.0X^) - a*Z^^(3a*x^.3X^)] , 
[(l-2v)a*2+l][^a*2F2(l)+F.,(a*)]+ ^ (a*2-i)[F.j(a*)-(l-2v) ^a*2F2(l)] 

(a*2-a *2)[(i-2v)a*2+i]. / (a*2-l)[(l-2v)a*2+a * 2 ] 

^ 0 
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CO-2v)a*2+a^*2][^ a*2F2(l)+F,(a*)]+(a*2-a^*2)(;p^(a*).(1.2v) ^ 

* r 

(a*2-a *2)[(i.2v)a*2+l> (a*2.i)[(i-2v)a*2+a *2] 

O 0 

(24 a-d) 


a*2F2(l)] 


and, for simplicity, it Is assumed that * V 2 « v. 

If one further assumes that = E 2 = E as well as * vg * > (i.e., 
that the elastic constants but not the thermal coefficients of the two 
cylinders are the same), then the expressions for hoop stresses are 
further simplified and become 


p 2 +a 2 a b 

I 0 -|(r,t)rdr + ag 02 (r,t)rdr] 

**n if a 


+ j e^(r,t)rdr - c*i,r2e^(r,t)} , 


(25) 


a. 


cTQ2(»'»t) ** ( b ' 2 -a° 2 - J ei(r,t)rdr + «2 J 

® a a 

0 


d r 

+ a-j I e.j(r,t)rdr + “2 | 


If the thickness a-a^ of the inner cylinder is indeed very small, 
then the analysis given in this section may be simplified quite consid- 
erably by using the identities 

Zoq(x,x) = 0, Z^ij(x,x) = 0, Z^i|(x,x) ^ • (27) 


In this case the characteristic equation (21) giving the roots may 
be approximated by 
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(28) 


An interesting special case of the thermal shock problem under 
consideration would be the yielding of the inner cylinder which may 
occur if the step change in the inner wall temperature is sufficiently 
high. The limiting value 0 p » which corresponds to the yielding 
of the membrane cladding may be obtained from (22) by letting 

cf0^(a*j+O) * cTp I (29) 

where op is a measure of the yield condition of the material (e.g., 
the flow stress). 

For the very thin inner cylinder mentioned above it can be shown 

that 

F,(a*) s 0 . 9,(a*,t*) = . 

and from (29) we find 
(l-v)cfp 


(30) 


(31) 


3. The^ Crack Problem 

The thermal stress problem for a composite cylinder containing an 
axial crack in a radial plane (Fig. 1) may now be solved by using the 
equal and opposite of the hoop stress given by (23) as the crack surface 
traction. As shown in [6], this problem must be solved under the 
following boundary conditions: 

cTj,i^(b » 0 »t ) “ Oj c7^g(b»0,t) — 0» (O^0^TT» 0<t) } (32ayb) 
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,t)-AagQ(a,e,t) - 0 , (0<e<27T, 0<t) , 

<jj,g(ai0,t) + X cfgg(a,e,t) * 0 , (0£e<2Tr, 0<t) 

* (e<r<g, 0<t) , 

U0(>^»O,t) = 0, (a<r<e, g<r<b, 0<t) . 

where cr^j. (1,j=r,e) and Ug are the stresses and the circumferential dis- 
placement in cylinder 2 for the perturbation problem and 

a 

"IT plane stress 

1 ,1 ®o> aE^d-vj^) 

' a^ aE2'(l-v,i!)+(ira^)e^v2(l+v2) ’ P’®"® strain. ( 36 ) 

By using the model to treat a cracked cylinder reinforced by a 
membrane developed In [6], the present transient thermal stress problem 
may be reduced to the following Integral equation: 

9 g 

1 + I k(r,p)f(p,t)dp. Ii^aj2(r.t) . 

e e 

(e<r<g, 0<t) , (37) 

where k = 3-4vg for plane strain, k = (3-v2)/(l+v2) for plane stress, 

y - E2/2(1 +v 2), and the unknown function f and the kernel k are defined 
by 

~ CUgCr.+O.t) - Ug(r, -O.t)], (e<r<g) , (38) 

k(r,p) = - 1 {- ^ + 2 yq + + 65^r 


(33) 

(34) 
(35a) 
(35b) 
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®o * i^TT-xTW^ XT (p ■ » 


^0 “ aT{i.^-y;.;bZ(y+xy 


C-(ln)p + . 


®i * rrraWn+xT ^ • 


. b-2(n-l) ^n-1 . ^ J1 , 

''0 0 

■^n * W C-a^*''^’V'*"'^’*(l-2Xn-x)-a2(l+n-nl-\) a^^"'^’'(l-2nx-x) ^] , 


c. = n 2 (l+ 2 x+A 2 )a’“b 2 - 2 (n 2 -l)(l-A 2 )-[l+ 2 nX+( 2 n- 1 )x 2 ]a’^'’b^" 
-r.V2n!.42fi-H)x2]a^V2'’ + n2(1-2x-3x2)a2b'2 , 
c„ = {(l-n2)(l-x2)-[l-2nx-(2n+l)x2]a^''b‘^'’+n2(l-2x-3x2)a2b‘2}p‘*"*’' 
+{(2+n-n2)(l-x2)b"^-(2+n-n2)(l+x)2a'^}p‘*"‘’' 
+'{(l+n)[l+2nx-(1-2n)x2]a‘^'’-[(l+n)-(n2+l)x23b'2'’}p''‘'' 

+ {(2+n)n+x)V^b'2"-(2+n)[l+2nx-(l-2n)x2]a'^'’b'2}p"*’ , 


-n- 


<!„ ' {-(1-n)[l-2nx-(1+2n)x2]a^'’+(l-n)(l->2)b^")p'*"*’^ 
.{(2-n)(l+x)2a‘V"-(2-n)n-2nX-(l+2n)x2]a^"b'^)p't""” 

■•■{n^O+x) W+0-n2)(l-x2)-n-2nx-(l+2n)\2]a^"b‘^")p'’‘’ 

+{(2-n-n2)n+x)2a'2-'2-n-n2)(l-x2)b‘=^)p'’+’ . (40) 

Fron (35b) and the definition of f as given by (38) it is seen 
that the integral equation (37) must be solved under the following 
single-valuedness condition: 

,9 

f(r,t)dr = 0 . (41) 

I 

e 

The singular integral equation may easily be solved by using the technique 
described, for example, in [7]. If the crack is fully embedded in 
cylinder 2 (i.e., if a<e<g<b), then the solution of (37) is of the form 

f(r,t) " , (e<r<g, o<t) (42) 

[(r-e)(g-r)]'^ 

where F is a bounded unknown function. In order to have a solution such 
as (42) it is, of course, assumed that the combination of mechanical and 
transient thermal hoop stresses in the neighborhood of the crack are pre- 
dominantly tensile. Otherwise the crack may close at one or both ends 
and (42) would not be valid. If part of the crack lies in a region of 
compressive hoop stress, then the crack surfaces may close smoothly and 
the location of the related crack tip would be unknown. In this case 
the unknown crack tip location is determined from the cusp condition 
which requires that the corresponding stress intensity factor be zero 
and the problem may be solved by using the technique described in [8]. 


After solving (37) the Mode I stress intensity factors at the crack 
tips e and g may be determined from 


k,(e.t) » ^ 11m f(r.t) » ^ 

' ' ^ r->e ' /(g-e)72 

k,(g.t) = - >^TFiT nr,t) . - • («) 

In the special case of the inner crack tip touching the reinforcing 
membrane (i.e., for e=a) it was shown analytically in [6] that at this 
crack tip the derivative f of the crack surface displacement is bounded. 
Physically this follows from the fact that since the membrane has no 
bending stiffness, when the crack touches it, it would form a kink. In 
this case the solution of (37) is of the form 

f(r,t) = , (a=e<r<g . 0<t) , (45) 

(g-r)'" 

where the bounded function F is determined by following the technique 
used in [6]. 

In the other special case in which the external loads are suffi- 
ciently high so that the membrane is fully yielded (if, for example, 
®ao^®F mechanical loading is present), the problem may be treated 

as an "edge crack" problem and the yielding membrane may be replaced by 
constant tensile tractions on the crack surface. In this case, too, the 
only relevant stress intensity factor is that at the crack tip r=g. 

Finally, there is the case of the broken membrane in which the 
problem may be treated as an edge crack problem. 

4. Results and Discussion 

Some calculated results giving the transient temperature and ther- 
mal stress distributions obtained from (16), (17) and (26) are shown in 


Figures 2-5. In (16) and (17) fifty terms are used in the series to cal- 
culate the temperatures which proved to be quite sufficient. The 
results are calculated by varying dimensionless quantities defined by 
(19). In these calculations it is assumed that for the two cylinders 
the elastic constants are equal (E^^Eg-E, ) but the thermal con- 

stants are different (D^?* 02 » k^ji^kg) and the composite cylinder is under 
a state of plane strain. Figures 2-5 snow the temperature and the hoop 
stress distribution for two different clad thicknesses ((a-aQ)/(b-a) * 
0.016 and 0.008). Referring to (15) note that for sudden cooling from 
inside is negative resulting in tensile stresses at and near the 
inner boundary (Figures 3 and 5). For an embedded crack (a<e<g<b) the 
corresponding stress intensity factors are shown in Tables 1 and 2. The 
tables show the normalized stress intensity factors which are defined by 


, k,(e,t) 

k,(e,t) = lim /2(e-r) a (r„0,t) , k(e) = -4 . 

' r-*e “o 

- ki(g.t) 

ki(g.t) = lim </2(r-g) a„.(r,0,t) , k(g) = — r , 

’ r-g ee 


Eoa>i0.^ 

^0 * - ''^ 9 ^ • 


(46) 

(47) 

(48) 


Note that for longer crack lengths the outer tip r=g of the crack would 
be in the compression zone and consequently the related stress intensity 
factor becomes negative. 

In the case of the inner crack tip terminating at cylinder-reinforcement 
interface the stress intensity factor at the outer tip r*g is shown in 
Figure 6. It may again be observed that as g increases the stress intensity 
factor k^(g,t) decreases and eventually becomes negative. In these cal- 
culations, too, it is assumed that E-|=E 2 =E, v-j=v 2 =v. 

The effect of modulus ratio E-|/E 2 =E* on the transient thermal 
stress cPg 2 in the main cylinder is shown in Figure 7 for the case of 
plane strain and in Figure 8 for plane stress. The figures show the 
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results at a fixed time t*=D-jt/b2=0.0005 and for 

Dg/D-j * 3 and kg/k^ * 3. Tables 3 and 4 show the corresponding normalized 
stress intensity factors for an embedded crack. The results for the 
crack touching the membrane and for the broken membrane (i.e, for the 
edge crack) are given in Table 5. Note that for the broken membrane 
case generally the stress intensity factors are considerably higher. 

Tables 1-4 and Figure 6 clearly show that an internal crack in a 
cylinder undergoing thermal shock at the inner boundary would always 
tend to propagate towards the inside wall, it may be arrested at the 
cylinder- reinforcement boundary if the reinforcement material is suffi- 
ciently tough, and under thermal shock alone it is not possible for the 
crack to propagate through the entire wall thickness of the cylinder. 

The tables also show that, generally for a given crack geometry the stress 
intensity factors first increase and, after going through a maximum, 
then start decreasing with increasing time. This behavior as well as the 
fact that for relatively long cracks the stress intensity factors at the 
two crack tips may have opposite signs is suggested by the thermal stress 
profiles shown in Figures 3, 5, 7 and 8. 
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Table 1. Normalized stress intensity factors for an embedded crack in a 
reinforced cylinder subjected to transient thermal stresses 
= 0.7968, I = 0.8, k2/ki -3,^=3,^= 0.8, | = 0.84, 

r .c .. ... 1 


Ei=Ej, 




SI 

t*«0.0005 

t*=0.001 

t*=0.005 

k(e) 

k(g) 

k(e) 

k(g) 

k(e) 

Mg) 

0.05 

0.11253 

0.07698 

0.16846 

0.13549 

0.11857 

0.10180 

0.10 

0.09617 

0.03173 

0.15356 

0.09109 

0.11120 

0,07845 

0.20 

0.06695 

-0.03653 

0.12573 

0.01508 

0.09707 

0.03509 

0.30 

0.04130 

-0.08115 

0,09912 

-0.04545 

0,08270 

-0.00419 

0,40 

0.01857 

-0.10986 

0.07337 

-0.09284 

0.06770 

-0.03938 

0.50 

-0.00189 



-0.13062 

0.04836 

-0.13134 

0.05204 

-0.07128 










Table 2. Normalized stress intensity factors for an embedded crack in a 
reinforced cylinder subjected to transient thermal stresses 
” 0.7984, = 0.8, 1^ ” 15” “ ^ §■ * 0,84, 


E, - E. 


= V,). 


g-e 

t *=0.0005 

o 

II 

■It 

*» 

001 

t **0.005 

b-a 

01 

k(g) 

k(e) 

k(g) 

k(e) 

k(g) 

0.05 

0.13677 

0.09793 

0.18877 

0.15414 

0.12371 

0.10688 

0.10 

0.11899 

0.04822 

0.17328 

0.10742 

0.11642 

0.08352 

0.20 

0.08715 

-0.02784 

0.14430 

0.02677 

0.10247 

0.04005 

0.30 

0.58931 

-0.07845 

0.11639 

-0.03816 

0.08822 

0.00057 

0.40 

0.3376 

-0.11135 

0.08916 

-0.08943 

0.07328 

-0.03486 

0.50 

0.01106 

-0.13495 

0.06259 

-0.13114 

0.05764 

-0.06692 



Table 3. The effect of the modulus ratio stress Intensity factors 

In a reinforced cylinder, » "h^“l^2®ao , a^/b * 0.7968, 

a/b « 0.8, e/b * 0.81, t* * Dt/b^^ « 0.0005, ag/a, - 0.8, kg/k^ « 3, 
Dg/D^ * 3, s vg = 0.3, the case of plane strain. 


gf! 

E,/E. 

j = 1 

£,/Ej 

! " 2 

E,/Eg * 6 


= 10 

kn(e) 

kn(9) 

k„(e) 

k„(g) 

k„(e) 

k„(g) 

“n<«> 

kn(g) 

0.05 

0.57556 

0.49535 

0.55469 

0.47486 

0.48101 

0.40217 

0.41844 

0.34024 

0.10 

0.55370 

0.39808 

0.53000 

0.37589 

0.45059 

0.29988 

0.38554 

0.23669 

0.15 

0.53438 

0.31039 

0.50762 

0.28680 

0.42198 

0.20803 

0.35420 

0.14377 

0.20 

0.51465 

0.23159 

0.48504 

0.20706 

0.39366 

0.12643 

0.32335 

0.06143 

0.25 

0.49342 

0.16180 

0.46139 

0.13673 

0.36513 

0.05505 

0.29263 

-0.01036 

0.30 

0.47052 

0.10099 

0.43654 

0.07569 

0.33636 

-0.00639 

0.26207 

-0.07195 



Table 4. The effect of the modulus ratio on the stress intensity factors 
in a reinforced cylinder, a^j/b * 0.7968, ‘ 

a/b = 0.8, e/b = 0.81, t* = Dt/b2 = 0.0005, ag/a^ « 0.8, kg/k^ “ 3, 
Dg/D^ = 3, = vg - 0,3, the case of plane stress. 


i 

E,/E. 

> = 1 

Ei/E. 

> " 2 

E,/E. 

> * 6 

E,/E2 

= 10 

■ 

kn(e) 

k„(9) 

kn(e) 

k„(g) 

k„(e) 

k„(g) 

1 

k„(e) 

k„(g) 

0.05 

0.40280 

0.34667 

0.38800 

0.33217 

0.33564 

0.28058 

0.29124 

0.23668 

0.10 

0.38732 

0.27847 

0.37022 

0.26256 

0.31276 

0.20793 

1 

0.26620 

0.16285 

0.15 

0.37362 

0.21700 

0.35405 

0.19996 

0.29130 

0.14295 

0.24263 

0.09703 

0.20 

0.35965 

0.16179 

0.33784 

0.14402 

0.27048 

0.08554 

0.22017 

0.03924 

0.25 

0.34468 

0.11292 

0.32100 

0.09474 

0.25002 

0.03562 

0.19862 

-0.01068 

0.30 

0.32858 

0.07035 

0.30344 

0.05204 

0.22987 

-0.00710 

0.17791 

-0.05312 


























Table 5. The effect of E-j/Eg on the normalized stress Intensity factor 
kj^Cg) for a crack touching the membrane cylinder inter'face 
and for a fully broken membrane, k^(g) « -kiCgtt)/E2a^9gQ/(g-a), 
a^/b « 0.7968, a/b=0.8, e/a=l, t*=D^t/b2=0.005, k?/ki=3, 

Dg/D^ * 3, = Vg = 0.3, 02/«i * the case of plane strain. 



Crack touching the membrane 

broken membrane- 
edge crack 

E,/E2=1 

E,/E2 = 2 

El/Eg - 10 

E,/E2 = 1 

E,/E2 . 2 

0.10 

0.23866 

0.22678 

0.15943 

0.39231 

0.39201 

0.20 

0.20237 

0.18697 

0.11076 

0.35107 

0.34859 

0.30 

0.16164 

0.14404 

0.06290 

0.32931 

0.32314 

0.40 

0.11804 

0.09940 

0.01651 

0.31389 

0.30323 

0.50 

0,07348 

0.05484 

-0.02742 

0.29972 

0.28361 

0,60 

0.02973 

0.01190 

-0.06839 

0.28100 

0.25851 

0.70 

-0.01222 

-0.02881 

-0.10741 

0.24839 

0.21943 
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Figure 3. 


Hoop stress in a hollow compos 
sudden change in the inner 
( (9"3Q)/(b-a) = 0.016 » ao/'^i ' 
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Transient temperature distribution in a hollow composite 
cylinder due to a sudden change in the inner surface tem- 
perature ((a-aQ)/(b-a) = 0 . 008 ). 
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Figure 5 . Hoop stress in a hollow composite cylinder caused by a 

sudden change in the inner surface temperature ((a-aQ)/(b-a)= 
O.OOSj 0^2/011 ” 0 . 8 ). 
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Figura 6. The normalized stress intensity factor k(g) = 

”5<(g,t)/Ea^0gQ /(g-e)/2 /(1-v) in a reinforced cylinder 
which contains a crack terminating at the reinforcement- 
cylinder interface and which is subjected to a sudden 
change in the inner surface temperature (E^=E 2 =E, v^=v 2 =v, 
(9”aQ)/(b-a) = 0.016, the case of plane strain). 
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Figure 7. 


The influence of the modulus ratio distribution 

of hoop stress in a composite cylinder undergoing a sudd#.n 
change in the inner surface temperature (t*=D^t/b^=0.0005, 
(a-ao)/(b-a) = 0.016, e/b = 0.81, the case of plane strain). 
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Figure 8. The influence of the modulus ratio E^/E 2 on the distribu- 
tion of hoop stress in a composite cylinder undergoing a 
sudden change in the inner surface temperature (t*= 
Dit/b2=0,0005, (a-ag)/(b-a) = 0.016, e/b = 0.81, the case 
of plane stress. 
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